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Abstract. We consider a free charged particle interacting with an electromagnetic 
bath at zero temperature. The dipole approximation is used to treat the bath 
wavelengths larger than the width of the particle wave packet. The effect of these 
wavelengths is described then by a linear Hamiltonian whose form is analogous to 
phcnomenological Hamiltonians previously adopted to describe the free particle-bath 
interaction. We study how the time dependence of decoherence evolution is related 
with initial particle-bath correlations. We show that decoherence is related to the 
time dependent dressing of the particle. Moreover because decoherence induced by the 
T = bath is very rapid, we make some considerations on the conditions under which 
interference may be experimentally observed. 



PACS numbers: 03.65.Yz, 03.70. +k, 12.20.Ds 
1. Introduction 



Decoherence, that is the destruction of coherent phase relations present in the elements 
of the density matrix describing the system, is induced by the interaction between 
the quantum system and the environment. It appears to be relevant in fields ranging 
from quantum measurement pQ to classical-quantum interface [2] , quantum information 
theory and computation [3] and cosmology [3]. Recently experimental evidence of 
environment induced decoherence has also been reported El HI El El CHI HI] • Several 
analysis of decoherence processes have been reported for the case of a particle, either 
free or in a potential, linearly coupled to the environment modelled as a bath of 
harmonic oscillators at temperature T O EH El EH EB UZj • Usually, initial conditions 
are adopted in which the system and the environment are decoupled, the interaction 
being effective after some initial time t = ^Hl E3 I2IH EI] It has been shown that, 
starting with factorizable initial conditions, which corresponds to the absence of initial 
correlations between the system and the environment, it is possible to separate the 
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decoherence into two characteristic parts. The first, related to the thermal properties of 
the bath, has a typical development time which in some models goes like T -1 ; the second, 
related to the zero point fluctuations of the oscillators of the bath, has a characteristic 
development time independent of temperature [2*2] . 

The time development of decoherence is also affected by the initial presence of 
correlations between the system and the environment O E3 123 I2H 123 123 I2ZI- 
However to choose the amount of correlations present in the initial conditions and to 
determine their influence on the time development may be a delicate problem mainly 
when interaction with the bath is always present as in the case of a charged particle 
interacting with the radiation field. However when the initial time is taken immediately 
after a fast dynamical evolution, like immediately after a collision, it appears to be 
reasonable to consider conditions where the particle and the bath are not in complete 
equilibrium [2B], the extreme cases being either of completely factorized or completely 
correlated initial conditions. Non factorized initial conditions have in fact been adopted 
for particles interacting with a thermal bath corresponding to the condition where a 
position measurement has been made on the particle once it is in equilibrium with 
the bath at temperature T [13 [23 EZ[- Other types of initial conditions have been 
also considered with the particle subject to a potential abruptly modified at t = or a 
combination of this and of the previous ones [21]. The presence of entanglement between 
a particle and a bath at zero temperature has been in particular explicitly taken into 
account [25] . 

In the context of quantum computing it has been shown, that in an ensemble of 
two level systems interacting with a reservoir of harmonic oscillators, decoherence among 
the two level systems develops because of the buildup of correlations between each two 
level system and the environment [3111 • A similar mechanism has been suggested to 
occur also in the case of a free charged particle interacting with the electromagnetic 
field vacuum. Here the decoherence among different momenta of the particle wave 
packet develops because of the buildup of correlations between each momentum and 
the associated transverse electromagnetic field structure that is responsible of the mass 
renormalization |31j . 

To treat the case of a non relativistic charged particle interacting, within the 
electric dipole approximation, with the electromagnetic field at temperature T, both 
the Hamiltonian approach ^3] and functional techniques [T3 ITE] have been used. 
In particular development of decoherence has been studied for charged particles by 
examining the time dependence of the interference among two coherent wave packets. 

Recently diffraction experiments have been performed where interference among 
different wave packets of the same particle has been observed also for rather long 
times [3 13 13 E] There are however indications that starting from uncoupled initial 
conditions decoherence develops also when the environment is at zero temperature and 
usually it may occur faster than the typical times of motion JTUI221- It appears thus of 
interest to examine how the development of vacuum induced decoherence depends from 
the amount of correlation initially present between the particle and the environment. 



Initial correlations effects on decoherence at zero temperature 



3 



Our model consists of a non relativistic free particle linearly interacting with a bath 
of harmonic oscillators at zero temperature. Our results will in particular be specialized 
to the case of a charged particle embedded in the electromagnetic field modes represented 
by a set of harmonic oscillators. This specialization is obtained from the general case 
by choosing the appropriate form of the coupling constants. 

To study the effects of initial condition on the time development of decoherence 
among the momentum components of the same wave packet, the behavior of the reduced 
density matrix elements shall be analyzed. 



2. Model 



Phenomenological Hamiltonians, where the interaction is described by a linear coupling, 
have previously been adopted in the study of decoherence of free particles interacting 
with the environment treated as a bath of harmonic oscillators [T2"| ITU ITH] . 

Here we consider a non relativistic free charged particle, initially moving at a 
velocity v, interacting with the transverse electromagnetic field. 

Taking the electromagnetic field as a set of normal modes, each characterized by 
a wave vector k and a polarization index j, the potential vector in the Coulomb gauge 
and with periodic boundary conditions taken on a volume V is 

Mr) = E f SLe-** + h,^) , (1) 



where the polarization vectors (j = 1, 2), r is the particle position operator and 

& k j and & k j are the annihilation and creation operators of the modes that satisfy the 
commutation rules [a^j, aJ k , .,] = dkk'^jj 1 - 

The non relativistic minimal coupling Hamiltonian has the form 



H 



eA(r 



P 



+ E f&kHj&kj , ( 2 ) 

k,j 



2m 

where p is the particle momentum operator, mo is the bare mass and e is the charge of 
the particle. 

In the case of a free particle the dipole approximation, typically adopted for 
bound charges, may be also adopted if the linear dimensions of the wave packet are 
small compared to the relevant wavelengths of the field. This can be applied to 
our Hamiltonian (J2J) introducing a cut off frequency Q such that the corresponding 
wavelength still allows the application of dipole approximation. Moreover in order to 
adopt this approximation without having problems related to the distance covered in 
time by the particle, taking advantage of Galilean invariance of the non relativistic 
Hamiltonian, we put ourselves in the reference system comoving with the particle so 
that its average initial velocity is zero [17 . All the quantities we have introduced 
are obviously relative to this reference system. Within the dipole approximation, the 
operator r can be replaced by a parameter r indicating the wave packet position at 
time t. In absence of interaction, in the comoving reference system, t*o is obviously 
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given by the initial position of the particle. Substituting this parameter in Eq. (J2J), the 
Hamiltonian becomes 

k = in~ + £ h ^l^ - ■ A ^ > ( 3 ) 

where we have neglected the quadratic interaction term. In fact its physical origin is 
linked to the particle average kinetic energy due to the vacuum fluctuations |32j, and 
it is usually very small compared to the linear term. On the other hand it can also 
be exactly eliminated by a canonical transformation of the Bogoliubov Tiablikov form 
[35] . This transformation has in fact been used to eliminate an analogous quadratic 
interaction term in the phenomenological Hamiltonian used to describe a free particle 
coupled to a dissipative environment j!2j . 

The advantage of using the dipole approximation and of having a linear interaction 
in our system manifests itself in the fact that the Hamiltonian can be treated exactly. 
Moreover the Hamiltonian of Eq. (J3J) is formally equivalent to a model Hamiltonian 
previously used, in the context of quantum computing, to study the decoherence of an 
ensemble of two level systems coupled to a reservoir of harmonic oscillators [SO] • This 
will permit to develop a physical analogy between these systems. 

It is clear that the wave packet describing the free particle is subject to spreading. 
This limits the range of time where dipole approximation may be used [14J. This 
spreading of the wave packet has been also taken in consideration with regard to eventual 
problems in the definition of decoherence [54] . 

Using Eqs. (jSJ) and (JIJ the Hamiltonian can be written as 

H = E £r *p + E MLAi + E *v 9lj + *ye <fc - r °) , (4) 

P k,j P,k,j 

where we have introduced for the momentum operator the notation J^pP&p, with 
&p = \p)(p\ the projection operator on the momentum eigenvalue p, while gfe ■ are 
the coupling coefficients given by 



9kj = ~P e k,j — \ 77—- ( 5 ) 

Eq. (0J) expresses the charge-field Hamiltonian in its unrenormalized form. To obtain 
the renormalized form we introduce the physical mass m as 

— = - 1 + — , (6) 
mo m \ m ) 

where 5m represents the mass variation due to the coupling with the bath. The 
renormalized form is then 

h = E £ * P + E + E * P 9h Uy-*- r ° + h^ k r A + e !~ * P 

P k,j p,k,j V P 

(7) 



The Hamiltonian under the form or (|7|1 can describe a variety of physical systems 
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with an appropriate choice of the coefficients g v k -. In fact a form analogous to it has 
been previously used, both in its unrenormalized [T3t fTB] and renormalized form [T2*lll3| . 
with the appropriate g k j, to treat a particle interacting linearly with a bath. 

In the following we shall keep unexplicit the expression of the coupling coefficients 
glp so that some results will be true not only for the charge- field interaction but in 
general for any particle-bath interaction with linear coupling. The explicit charge-field 
form (jSJ) will be only used at the end to evidence the dependence of the results from the 
parameters of the system. 



3. Fully correlated states 

To treat the effects of the presence of particle-field correlations in the initial state, we 
consider at first an initial condition with the particle in complete equilibrium with the 
zero point field fluctuations and localized at the position r$. This implies that each 
momentum component of the wave packet is in equilibrium with the fluctuations of the 
electromagnetic field. The momentum of the charge commutes with the Hamiltonian 
(J7J), thus the equilibrium state of a given momentum must be an eigenstate of H. In 
order to find these eigenstates we diagonalize exactly H by a canonical transformation 
using the unitary operators 

DL = exp [E %^ (4,e-* fc ^ - h,^)} , (8) 



k.j 



that act in the particle-field mode Hilbert space TCs <8> Tip 
The action of the canonical transformation, induced by D\a on the operators &k,j and 
& kj , is 

p 

sr, . _ nP- 1 ^, .nP _ s , y (T P g kJ ik-r 



/' 



a 



4j - dIj - 4j + Yl ^J^ 3 e tk r ° ■ ( 9 ) 

To treat properly the initial condition of particle and bath when correlations are present 
at all frequencies, we apply the canonical transformation given in Eq. (jHJ) to all the 
modes. This amounts to express H in terms of the transformed operators &k,j and a£ -. 
Thus, inverting Eq. (jUJ) and substituting into the Hamiltonian of Eq. (J7J). we obtain 

P k,j P P,k,j ™ 

In the last term only er p is an operator thus it is possible to write 

y 9kj = P 2 Sm 
hcuk 2m m ' 

k,3 

by an appropriate definition of 5m. In particular for the charge-electromagnetic field 
case, using Eq. (JSJ for g k p taking the continuum limit, the mass variation 5m has the 
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form 

AahTLm 2 

6TXC Z ITIq 

where Q is an upper frequency cut off. We see that, within second order perturbation 
theory in the charge, it coincides with the usual electromagnetic mass variation due to 
the interaction with the electromagnetic field [35]. 
In this way, inserting Eq. (jllj) in Eq. (jlU|) . the Hamiltonian reduces to 

H = E ^ & p + E ^* s L s fcj • ( 13 ) 

The eigenstates of belong to the complete Hilbert space Tis ® 7~Cf and have the form 
\p, {nk,j}) — I\k,j ^k,j\p)\ n k,j)- 111 the case where each momentum is in equilibrium with 
the dressed vacuum state, we consider the states 



In Eq. (JT3j) each operator Z)[ p ] acts only on the field- mode Hilbert space Hp-' . Because 
of the form of D^j, the state \{otk,j(p)}) indicates the product of the coherent states of 
all the modes of the field, each of amplitude ockjip) — 9k,j exp(— ifc • r )/hu>k, associated 
to the component \p) of the wave packet in the space TCs- The state \p) represents 
in the electromagnetic case the state of the coupled system formed by the particle of 
momentum p plus the transverse photons associated to it. 
The action of Hamiltonian (jlHj) on the states \p) of Eq. ()14j) reduces to 

H\P) = ^\P)- (15) 
Thus \p) are eigenstates of H with eigenvalues p 2 /2m. 

In the Schrodinger picture the time evolution operator expressed in terms of the 
transformed operators is 



\p, fe}) = lPL-b>IM = \p) n^Sl°w> = \p)\Mp)}) ■ (14) 

k,j k,j 



U(t) = exp 



P k,j 



(16) 



We consider an initial localized wave packet whose components are the correlated states 
of the particle in equilibrium with the bath in its vacuum state 

l*> = (17) 

p 

The initial density matrix of the total system is then given by 

p(o) = Y,c p c;,\p)(p>\ = Y,c p c;,\p)(p>\ ® \{o k , 3 })({o k >, r }\ . (is) 
p,p' p,p' 

From Eq. (|18|) the reduced density matrix of the particle at time t is given by 

p s (t) =tr F {U(t)Y,C p C;,\p)(p'\®\{O k ^ (19) 

p,p' 
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Inserting the time evolution operator (J 16)) and the expression for \p) (|14j) in Eq. (|19|). 
the matrix elements of the reduced density matrix can be cast in the form 



fff (t) = C P C;, exp 



it (p 2 — p' 



r2\ 



tr F {\{a kd (p)})({a k , d/ (p')}\} . 



h 2m 

The time dependence of pg P (t) is given simply by exp [— it(p 2 — p' 2 )/2mh], and 
represents the free evolution of the initial reduced matrix elements 

f?/(0) = C p C;,tT F {\{a kj (p)}}({a kl Ap'W • (20) 
It follows that for an initial state consisting of a coherent wave packet of particle-bath 
correlated states \p), in the form of Eq. (|THj) . decoherence doesn't depend on time. The 
decoherence present in the reduced density matrix elements is contained in the factor 
trp {\{cv k} j(p)}){{a k iji(p')}\} appearing in Eq. (j^UJ). For charge-field interaction, this 
factor can be interpreted as a consequence of the cloud of transverse virtual photons 
being associated to each momentum component of the particle. 

This interpretation of the presence of decoherence in our system is analogous to the 
one given for the case of two level systems linearly coupled to a bath of harmonic 
oscillators [SO] , where decoherence is attributed to the buildup of correlations between 
the environment and the states of the two level systems. 

By performing explicitly the trace on the field in the factor appearing in Eq. ()20|) 
we obtain 



tr F {\{a kJ (p)})({a kl , f (p')}\} 



k",j" n k,k',j,j' 



Y[ (0 fc /y|exp 

k,k',j,j' 



P 

9 k > j' 



huj k > 



k'-r 



k',j 



/C 



ik'-r 



x exp 



n ( 



p 

9l 



tlUJ k 
p' p 

9 k >,j> i 9 k j 
Tiujui hui 



t e 



ik-r s jk-r 

— d, k jV 



|0*j> 



>■ 



(21) 



k,k',j,j' '"^k' iiMJ k 

Hence, the trace reduces to a product of coherent states. We calculate Eq. (|2"T]) explicitly 
in the case of a charged particle interacting with the electromagnetic field. We use the 
expression for the of Eq. (JHJ), exploit the form of the scalar product between coherent 



states and perform the continuum limit of the field modes, J2h 



(2ttc) : 



Jo 00 d 3 u;, to get 



teF{\{a k j(p)})({a k ,, f (p')}\} 

, 2" 



Il eX P 

exp 



p p 
9 k ,j ~ 9 k ,j 

2h 2 u 2 k 



J|exp 

k,j 



e 2 2nh [{p-p') -Skj] 



[ mg Vu k 
P') " £ k,j] 2 



2h 2 u; 2 k 



(22) 



/|2 



exp 



a\p — p 

37imQC 2 



du;- 



-ui/Si 



37r 2 c 3 mo^ Jo uj^ 

where a = e 2 /hc. The last equation has been obtained summing over the polarizations 
and introducing as usual a high frequency cut off factor in the integral over frequencies 
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[To] . However the frequency integral appearing in the last equation maintains the 
infrared divergence. This makes the reduced density matrix of the particle diagonal 
and the decoherence between different momentum states of the particle is therefore 
complete. Thus, as a consequence of the infrared divergence, there is a super selection 
rule in the momenta with the emergence of stable super selection sectors [HE]- 
However, by taking into account the finite time measurement of the particle momenta, a 
lower cut off frequency w may be introduced, the frequency w representing the resolution 
in the detection process. Introducing this lower frequency cut off w in Eq. (j22J) we obtain 



{\{ak,j(p)})({a k i,ji(p')}\} = exp 



exp 



a\p — p 

3nmlc 2 
a\p — p'\ 2 

37rmlc 2 



du- 



CO 



(23) 



7 + In 



E 

n=l 



n n 



l) n (w 



where 7 = 0.577216 is the Euler's constant. 
For vj/VL <C 1 Eq. ("""«")) can be approximated as 



./12 



n r / m \ / r / m n \ a \P ~ P i w 

{ I {a k ,j (P) }) {{aw (P ) } I } = exp j 3nm 2 c 2 ln 7^" 



and the initial reduced density matrix elements of Eq. ()20|) take the form 

1 1 ■ \p - i> ' 1 ' 



pT'(0)^Wexp|^-|Xln^ 



(24) 



(25) 



Now decoherence between momenta in the matrix elements of Eq. ()25j) is not complete 
anymore and depends both on the distance between the values of the particle momenta 
and on w. 



4. Partially correlated states 

The need of introducing a lower cut off frequency can also derive from the initial packet 
preparation rather than from the final measurement. If the initial packet is localized 
by a measurement with an uncertainty Ax and then with a momentum uncertainty 
Ap ~ h/Ax, this is equivalent to a minimum measurement time At ~ h/ApAv [2*2*1 138j. 
to which it is possible to associate a lower frequency limit w ~ 1 / At to modes that have 
reached equilibrium with the particle. 

This indicates that we can start from an initial condition with the particle in equilibrium 
only with the modes at a frequency uj > w while it has no correlation with those at 
uj < w. Thus we take as initial condition a state whose momentum components have 
the form 

ro/c 00 

ip-> = b)nniM n hi)- ( 26 ) 

j k=0 k=zu/c 
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From Eq. we build the initial wave packet ^2 p C p \p m ) that corresponds to the initial 
density matrix 



vj/c 



(27) 



^(o)=Ewip}( P 'inniy^i n fe>fei 

p,p' 3 k=0 k=zu/c 

In Eq. f)27j) the modes at frequency higher than w are distinguished from those at lower 
frequency and we have eliminated the cross terms between the two mode regions. In 
fact these terms won't contribute to the reduced density matrix elements because they 
will be eliminated in the field trace and then we don't consider them. 

To discuss the time evolution from the initial state represented by the density 
matrix p ro (0) it is useful to separate in the Hamiltonian of Eq. ((7J) the parts at frequency 
larger and smaller than w, i.e. 

J p 2 5m 



+ y <r p 

„ 2m m 



y 2m °" p 



zu/c zu/c . 

j,k=0 P,j,k=0 V 



(28) 



+ 



Y huj k &i ;j & ktj + Y a ? dij \ H,j e 

j,k=zo/c p,j,k=zu /c 



Ep 5m A 
<r„ . 
„ 2m m 



To treat properly the presence in the initial condition of correlations with the high 
frequency modes, in H, we apply the canonical transformation (jHJ) only to the modes 
at frequency higher than w. This amounts to express only in terms of the 

transformed operators H k j and a| -. Following the same steps leading from Eq. ((7|) 
to Eq. Jin| we get 

p2 oo ^2 fi m oo gP2 

H = Y 7T~ & p + ^< ro + Y h^k&kj&kj + Yo & p - Y T^Vv ■ ( 29 ) 

P j,k=w/c P p,j,k=m/c 

In the last term of Eq. (f2"TJj) only a p is an operator and, in analogy to Eq. (JTTJ) , it is 
possible to write by an appropriate definition of 5m >VJ 



9 k j _ P 2 Sm >z 



Y 

j,k=-€u/c 



%UJ k 



2m m 



(30) 



Sm >zu corresponds to the mass variation of the particle due to dressing with the high 
frequency modes. Thus, for the sum of the last two terms present in Eq. ()29j) we obtain 



j^5m A ~ 
^2m m ap ^ 



9k j P 2 $ m <™ ■ 
(Tp = 2^- a 



2m m 



v ■ 



(31) 



P "" p,j,k=-cu/c '"™ K P 

where 5m < - DJ = 5m — 5m >zu indicates the mass variation of the particle due to dressing 
with the low frequency modes. Using Eq. (|31|) in Eq. (|29|). the Hamiltonian H becomes 



H = Y 7T~ °P + Y o ~ °P + Hu<vj + H U>VJ = Y 

p 2m p 2m m p 



7T— + H t 

2m 7T7 



(32) 
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where H, 




and the mass m TO is defined by 



(33) 

No renormalization term is present in the Hamiltonian of Eq. fl32(l . This is due to 
the fact that only the interaction with modes at frequency u < w appears in H and 
that the mass m ro represents the bare mass with respect to the low frequency modes. 
The Hamiltonian being now expressed in terms of the bare mass, the low frequency 
renormalization term disappears. 

In the Hamiltonian H (p?2J) the terms S p (p 2 /2m ro ) & p , H^ <m and H w>m commute 
among them. This allows to analyze in a transparent way the influence of each of 
these terms on the time evolution of the reduced density matrix elements of p m (0) 
and to discuss separately the effects of high frequency correlated and low frequency 
uncorrelated modes. In fact we can write the reduced density matrix elements of the 
particle at time t as 

fflZ® = (p\trAUUt)pUW^(t)}\p') , (34) 

where U w (t) is the time evolution operator, that may be written as the 
product of three operators U m (t) = Us,m{t)Uu}<za{t)U u>w {t). The first, U s ,m = 
exp —iJ2 P (p 2 /2m m h)a- p t , depends on the kinetic energy and can be shown to affect 
only the free evolution of reduced density matrix elements. Introducing the time- 
ordering operator T<_, the second term is U w<m (t) = T^exp — % J * dsiT w<ro (s)//?, , and 

depends on the modes at frequency lower than w, while the third term is U u>ZJ (t) = 
T^exp —i J dsH UJ> ^(s)/h and depends on those at frequency higher than w. Thus, 
Eq. (134)1 can be written as 

it (p 2 — p' 2 ) 



P P sl(t) 



exp 



h 2m T 



ip\tr F {U^{t) U^(t)pU0)U-'(t) U-'(t)}\p') -(35) 



U w<w acts only on the low frequency part of the density operator p ro (0) of Eq. (|27j) 
representing the modes uncorrelated to the particle momenta. This leads to a time 
evolution of this part of p ro (0) formally identical to the one we would get from a fully 
uncorrelated initial condition. 

For a charge interacting with the electromagnetic field we may directly use, for the 
evolution of the low frequency part of the results already obtained in j^I], that 



is: 



P7 (t) = &* (0) exp 



it (p — p ) 



expT p ' p '(t) 



(36) 



h 2m 

where the complex decoherence function T p,p '(t) = T p,p ' (t)+iT p ' p (t) has been introduced 
with TP'P(t) and T p ' p (t) given respectively by the integral expressions 



r p /{t) 



2a\p — p 



t\2 



3tt ttiqC 



2^2 



da; exp 



(1 — COS Ult) 



(37) 
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which corresponds to the expression of T p ' p '(t) given in [SI] in the limit T = 0, in which 
only the vacuum contribution remains so that T p ' p ' (t) = r^ c '(t), and 

Here, because of Eq. (l2T|) . in the integrals defining T p,p (t) and Y p,p (t) we replace the 
cut off factor exp (—u/Q) with the step function Q(w — uj). In such a way we obtain a 
new expression for the real and imaginary part of the decoherence function, T p,p ' c {t) and 
T p ' p (t), which will depend only on the modes at frequency lower than w as: 

ppjZ/.n 2a \p-pf [°° (1 - cos ut) 



2a\p — p 



37r mpC 2 



/|2 

[-/ -Ci{wt) + In wt} , (39) 



ppj//.n 2e 2 (p 2 - p' 2 ) /•« / sinwA 2a p 2 -p' 2 



37r/?,moC : 

where Ci(^i7t) and Si(wt) are the cosine and sine integral function [3*7] . 

Using the expression of Ci(wt) for small and large values of wt [37] we obtain for T p,p ,{t) 

2a \p-p'\ 2 w 2 t 2 



37T TUqC 2 4 

2a\p — p 



3n m\c 2 



/|2 

In wt for wt ^> 1 



for rot 1 

(41) 



Eq. f)41|) shows that there are a quadratic and a logarithmic time evolution regimes. 
The transition from the first to the second one occurs at a typical time t pa w~ x with 
|r^'(f)| = Q ~ |f ^-j^ ■ ■ The effective magnitude of T p ' p ' c (t) at a given time depends 
both on w, which is determined by the preparation or observation procedure, and on 
Q which depends on the charge and the mass of the particle and on the distance of 
the reduced density matrix element from the diagonal in the momentum space. The 
vacuum influence on the time evolution of decoherence for Q <C 1 remains small because 
of the logarithmic dependence of Y P)P ' c (t) from time. Its effect increases with Q and for 
sufficiently large Q decoherence may become effective at times t > w^ 1 . This behavior 
is shown in figure 1. 

In conclusion we see that the effect of the low frequency modes, described by U^^t), 
reduces to a multiplication factor in the reduced density matrix element p^^if) of the 
form exp [f p £(t) + iT p ' p ' (t)] where T p £{t) is given by Eq. flU) and V p ' p \t) by Eq. flU). 

With regards to the third term in the Hamiltonian (|32jl . H U>ZJ , representing the 
high frequency modes contribution, its action on the time evolution of p^(t) can be 
obtained by following the equations that from Eq. (jlfijl lead to Eq. (j25|) . where the lower 
frequency represents in that case the effect of the finite resolution of the detection 
process. We are therefore led to a term formally identical with the one of Eq. (|25jl that 
doesn't influence the low frequency part of p ro (0). 



Initial correlations effects on decoherence at zero temperature 



12 




Joining together all the previous considerations we finally obtain for the complete 
expression of the reduced density matrix elements 



p p sl(t) = C P C;, exp 



x exp 



it (p 2 — p' 2 ) 



k 

2 a \p 



2m T 
P 



exp 



.2a (zut — Si(wt)) 

StIITLqC 2 



(P 2 -P' 2 ) 



'12 



37rmoC 2 



[7 — Ci(zut) + In tut] > exp 



a\p — p 

'i r nm\c 2 



In 



(42) 



We observe that with regard to the contribution in Eq. (|42jl of the low frequency modes, 
we have used Eqs. (|39|) and (|40|1 whose form is valid at all the times. For the contribution 
of the high frequency modes we have used Eq. (|25j) which is valid in the limit zu/Q <C 1. 
However it is possible from Eq. (|42|) to obtain a simplified form of p P g P m {t) for small and 
large tut, using for F^f c '(t) the corresponding expressions of Eq. (jJTJ). 

The form given by Eq. (|42jl describes the time evolution of the reduced density 
matrix elements when, in the initial state, correlations with the field modes at frequency 
larger than w are present. The form for ps^(t) of Eq. (|42|) is intermediate between the 
two extreme forms, completely correlated (|20 [ 122 (1 and uncorrelated (|36|h and it coincides 
with them respectively in the limits w — > and w — > fi. 

One must observe that in the fully correlated case {w = 0) the particle is completely 
entangled with the field and the initial reduced density matrix (jl8)l contains already all 
the decoherence, which therefore remains constant with time. For w 7^ the initial 
entanglement is partial and it is the coupling dynamics that induces the progressive 
entanglement with the modes at frequency less than w thus leading to an increase of 
the decoherence present at t — 0. 



5. Summary and Conclusions 



We have studied the decoherence among the momentum components of a free particle 
wave packet linearly interacting with a zero temperature bath by working out the time 
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dependence of the off diagonal elements of the reduced density matrix of the particle. 
We have also examined the influence on the decoherence development by the initial 
conditions and in particular by the presence of correlations between the particle and the 
environment. 

The influence of initial conditions on decoherence has been previously studied by 
examining the attenuation of coherence between a pair of coherent wave packets in the 
case of a particle interacting with a bath either at finite [27] or zero temperature ■ 

We have considered initial conditions ranging from the case of absence of initial 
correlations to the case where they are fully developed. The particle-bath system is 
described in the former case by a factorized density matrix and in the latter, which 
corresponds to a stationary condition, by a linear combination of the eigenstates of the 
total Hamiltonian. The reason to consider this range of initial conditions is linked to the 
fact that when the interaction between the system and the bath is always present, as in 
the case of a charged particle interacting with the radiation field, it isn't appropriate to 
start from a condition where no correlations are present. On the other hand because of 
either a finite measurement time or preparation time it is not appropriate to start from 
a condition where there is complete equilibrium among the particle momenta and all the 
field modes. By taking for example a finite preparation time r, one is led to consider a 
situation where equilibrium correlations have been established only with the modes at 
frequency higher than w ~ 1/r, while modes of lower frequency are uncorrelated. 

This intermediate initial situation may be described by a density matrix formed 
by two parts: the first, time independent, describes the initial correlations already 
established with modes at frequency higher than w; the second time dependent, will 
describe the establishing of correlations with the modes at frequency lower than w. We 
have shown that, for the case of linear interaction, the increase of decoherence is at 
first rapid and goes like t 2 while, after a transition time of the order of slows and 
goes as lnt, in a certain sense reaching a plateau at t ~ w" 1 . The effect of partially 
correlated initial conditions does show only in the transition time, while the value that 
decoherence reaches at the transition time is only a function of the physical parameters 
of the particle like the value of the coupling constant, the particle mass and the distance 
of the matrix element from the diagonal. 

We have explicitly obtained our results for the case of a charged particle interacting 
with the electromagnetic field. In fact, by an appropriate choice of the width of the 
particle wave packet we have shown that it is possible to adopt the dipole approximation. 
In this context and by considering the case of low field intensities, we have shown 
that the interaction Hamiltonian becomes linear. In this case then the value of the 
decoherence at the transition time w~ x depends from the combination of parameters 
given by Q ~ |a \P~P 2 1 m Moreover the time independent part of the reduced density 
matrix p ro (0) given by Eq. ()27j) presents an infrared divergence when the preparation 
time is large and thus w — > (^0^) . In this case the decoherence in momentum space is 
complete and the divergence is connected with the establishing of a complete correlation 
with the low frequency modes. For finite w, with the buildup in time of correlations 
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with the low frequency modes these become populated. This process gives rise to a cloud 
of photons around the particle which are then responsible of mass renormalization. The 
two phenomena of decoherence due to zero point modes and dressing of the particle 
are thus related. Moreover when one performs an experiment that takes a finite time, 
that shows the presence of decoherence, like in interference experiments |H1E], it is 
appropriate to ask the form of the reduced density matrix that correctly describes the 
experiment. From what said previously it appears that the trace must be performed 
only on those field modes where supposedly equilibrium hasn't been reached. Then the 
characteristic time in which decoherence becomes complete is a factor Q/w larger than 
the characteristic decoherence time one obtains starting from completely uncorrelated 
initial condition 31 . 

Our analysis has been conducted in the context of non relativistic QED which is 
an effective low energy theory with the cut off frequency Q, in the spirit of modern 
quantum field theory, parameterizing the physics due to the higher frequencies jSH] • For 
this reason our final results must show a dependence on Q, that is however as usual weak 
(logarithmic), as for example in the case of non relativistic expression for the Lamb shift. 
The appearance of the bare mass m in our results, e.g. in Eq. (J20|) and (jI2J), is due to 
its use in the definition of the coupling coefficients g k j (JHJ). The g v k ^ in Eq. ((7J) and (|28j) 
could be expressed in terms of the fully and partially dressed mass by introducing a 
further renormalization term which is however 0(e 3 ) and therefore typically neglected. 
This procedure would eventually lead to final reduced density matrix expressions given 
in terms of the dressed masses with a correction in the exponents of 0(e 4 ) with respect 
to the actually results. This may be alternatively seen by expressing, directly in the 
final expressions for decoherence, the bare mass in terms of the dressed masses obtaining 
a further term containing the product of the mass variation of 0(e 2 ) for a factor a. 

As a final consideration we observe that a state describing an electron not in 
complete equilibrium with its surrounding field was suggested by Feinberg [10] to 
represent the electron immediately after a scattering event of duration r. In fact, due 
to causality requirements the electron can reach equilibrium only with its field within a 
region of size I < cr. Similar states of incomplete equilibrium have been also studied in 
QED when rapid changes occur in atomic or molecular sources j2H]. However the most 
promising field where states of not complete equilibrium may be implemented is solid 
state physics. The use of ultra fast spectroscopic techniques has in fact permitted the 
generation of almost bare electron-hole states and to follow the time evolution of the 
dressing process due to the interaction with the phonon field Taking into account 
the previous considerations, one can envisage a process where, starting from states of not 
complete equilibrium, the evolution of decoherence gives rise, in principle, to observable 
effects. Let's in fact consider the scattering of a charged particle from two scattering 
centers, the process lasting a finite time r « za~ x . From the considerations developed 
in this paper one must expect that, after the scattering, the state of the particle may 
be described by Eq. (J27|) and the decoherence which develops after the scattering gives 
rise to a decrease of the interference effects which depends on w. 
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